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We apply the point form of relativistic quantum mechanics to develop a Poincare invariant 
coupled-channel formalism for two-particle systems interacting via one-particle exchange. This 
approach takes the exchange particle explicitly into account and leads to a generalized eigenvalue 
equation for the Bakamjian- Thomas type mass operator of the system. The coupling of the ex- 
change particle is derived from quantum field theory. As an illustrative example we consider vector 
mesons within the chiral constituent quark model in which the hyperfine interaction between the 
confined quark-antiquark pair is generated by Goldstone-boson exchange. We study the effect of 
retardation in the Goldstone-boson exchange by comparing with the commonly used instantaneous 
approximation. As a nice physical feature we find that the problem of a too large p-uj splitting can 
nearly be avoided by taking the dynamics of the exchange meson explicitly into account. 

PACS numbers: 12.39.Ki,21.45+v 
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I. INTRODUCTION 

In 1949 Dirac formulated a way of incorporating rela- 
tivity into quantum theory that differed from quantum 
field theory Q. Although Dirac's paper was written in 
the context of classical mechanics, his methodology - the 
use of representations of the symmetry group of the the- 
ory of special relativity, the Poincare group - was also ap- 
plicable to quantum theory and to quantum field theory 
(for a review, see Ref. Q). In Dirac's original presenta- 
tion he made evident how to add interactions to a theory 
of free particles in agreement with the Poincare algebra, 
ending up with conditions for the interaction terms that 
are in general nonlinear. In 1953, Bakamjian and Thomas 
gave a prescription for an explicit construction involving 
only linear constraints [3j. 

Of the various forms that Dirac introduced the instant 
form is the most widely used 0, IE IE 0j although al- 
most exclusively in the context of quantum field the- 
ory. The front form of Hamiltonian dynamics became 
popular as a natural framework for treating parton phe- 
nomena. For a topical review, see Ref. |8|]. The point 
form of relativistic dynamics has also been considered 
in quantum field theory [E flfl ITlj . but because of its 
complicated quantization surface it was not further de- 
veloped. Only recently has the point form been rediscov- 
ered, this time in the context of the quantum mechanics 
for finite degree-of-freedom systems. Lev has analyzed 
electromagnetic current operators in the point form and 
then transformed them to the other two forms E2 • Klink 
has suggested a basis of states, called velocity states [l3| . 
which are suitable for few-body point-form quantum the- 
ory. These states have also been introduced by Kar- 
manov in a different context [l4l ITU . Klink made use 
of velocity states when treating nuclear physics problems 
EE 13 El Recently, also the Graz group employed the 
point form of relativistic dynamics to describe the elec- 



troweak structure of baryons within a chiral constituent 
quark model El EE El. 

This paper uses the point form to elaborate on a 
coupled-channel formalism which is applicable to a wide 
range of problems. As a first and simple application we 
have chosen a two-particle system of one (constituent) 
quark and one (constituent) antiquark which form vec- 
tor mesons. The hyperfine interaction in this simple sys- 
tem comes from the chiral constituent quark model with 
pseudoscalar meson exchange. Such an interaction has 
been used in a semirelativistic form wit h g reat success 
for the calculation of baryon spectra |5E |23[ . For vector 
mesons [24| the results within this semirelativistic ap- 
proach are not as good. In the present paper a fully 
relativistic calculation with the exchange-meson chan- 
nel explicitly included is presented and compared to the 
semirelativistic approach in which the meson-exchange is 
treated in an instantaneous approximation. 

The necessary formalism is introduced in detail in 
Sees. HTHVI In Sec. [n] we summarize relevant features 
of the Poincare group. Velocity states are introduced as 
a suitable basis for the quantum-mechanical treatment 
of few-particle systems in point-form. These states are 
subsequently used to construct the elementary meson- 
(anti)quark vertex, which enters the invariant-mass op- 
erator. The mass operator is treated in Sec. 1 1 1 1 1 in our 
example its interacting part arises from a pseudoscalar 
Hamiltonian density. Our mass operator is of Bakamjian- 
Thomas type and acts on a Hilbert space which is the di- 
rect sum of two-particle and two-plus-one particle Hilbert 
spaces. Such an ansatz deals with effective degrees of 
freedom in contrast to a quantum field theory. It leads 
to the coupled two-channel problem outlined in Sec. IIVI 
The validity of Poincare invariance for systems with a fi- 
nite (but not necessarily conserved) number of particles 
is guaranteed by the Bakamjian-Thomas construction as 
described in Sec. IIVI The fact that quarks and anti- 
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quarks are always confined is accounted for by adding 
harmonic oscillator confinement terms to the square of 
the kinetic terms of the coupled-channel mass operator. 
As a first step to solve the eigenvalue problem for the 
mass operator the two-channel problem is reduced to a 
one-channel problem with an optical potential which de- 
pends on the mass eigenvalue to be determined. The 
harmonic-oscillator eigenfunctions of the pure confine- 
ment problem are then used as a basis for expanding 
the quark-antiquark wave functions of the full problem 
including the hyperfine interaction. As a result, one can 
discretize the dynamical two-particle equation and obtain 
a set of coupled algebraic equations that can be solved. 
The eigenvalues of the mass operator are determined by a 
resonance condition, which takes into account the nonlin- 
ear appearance of eigenvalues in the eigenvalue equation. 
The structure of the equation and the implementation of 
confinement are discussed in Sec. 

Sec. I VII contains some remarks on the numerics and 
the specific solution method employed in the calculation. 
Comments on the instantaneous approximation and the 
choice of the model parameters are given in Sees. IVIII 
and IVIIII respectively. The results of the calculations 
are vector-meson masses and some branching ratios of 
their hadronic decay widths. The numbers are presented 
in Sec. IIXI they are compared with experimental num- 
bers as well as with the instantaneous approximation to 
elucidate retardation effects coming from the hypcrfinc- 
interaction. Concluding remarks can be found in Sec.lXI 
Our conventions, normalizations and matrix elements re- 
quired for the calculations are summarized in the ap- 
pendix. 

II. POINCARE GROUP 

The starting point for dealing with few-body systems 
in relativistic quantum mechanics is the set of commuta- 
tion relations of the Poincare generators 

[PM = o, (i) 

[Jliu,P K ] = iiSunP^ - Q^Pv) , (2) 

gvxJ^n — 9)1\Jvk) ■ (3) 

One can write these relations in a global way by defin- 
ing U A as the unitary operator representing the Lorentz 
transformation A on the Hilbert space. In the point form 
all interactions are contained in the four-momentum op- 
erator, so the significant commutation relation is 

[P*,P1=0, (4) 

which states that the components of the four-momentum 
commute among each other. The other commutation re- 
lation involving P M is written as 

UaP^ 1 = (A-Yr , (5) 
which means that the four-momentum operator has to 
transform as a four-vector under Lorentz transforma- 



tions. The commutation relations of the Lorentz gen- 
erators among themselves are unaffected by interactions 
in the point form. We will refer to Eqs. (@J and © later, 
when we construct mass operators. 

Let us start with single-particle states \p, a), e.g. for 
spin-i particles, which transfer irreducibly under the 
Poincare group. Defining the action of the four-momen- 
tum operator P^ on such single-particle states by 

P' t \p,a) = \p,a)pf t (6) 

one can easily show that Eqs. and (JjjJ are satisfied 
for a single-particle representation. In the following we 
need also the Poincare transformation properties of such 
states: 

U b \p,a) = e- ibp \p,a) and (7) 
U A \p,a) = Yl \Ap,° f )Di,„(R w (p,A)) . (8) 

<x'=±± 

Ub denotes a space-time translation by a constant four- 
vector b and the Wigner rotation R\y (p, A) is given by 

R W ( P ,A) = B-\A^)AB(^). (9) 
m m 

i 

The D£, are the matrix elements of the standard Wigner 
IP-functions |25j . 

With Eqs. il(ili)l) we see that Eq. |@J| is immediately 
satisfied. For Eq. (J5J we have, applying its left-hand side 
to \p,a), 

U K P»Ul x \p,a) = 
= U A P»U A -i\p,a) 

= UkP^^P^DI'ARw^A- 1 )) 

a' 

= U A J2 |A- X P, a'^A-'r^D^Rwip, A" 1 )) 

= U h U A -Ap,a){A- l y u p v 
= U^-^a^A-y^ 

= (A-y„P>,a) , (10) 

which is the desired expression for the right-hand side. 
In this derivation we have extensively used the represen- 
tation properties of U A and U^ 1 , e. g. that 

UaUa-i = U AA -i = Ut = l. (11) 

In the generalization of single-particle states to mul- 
tiparticle states, it is useful to introduce velocity sta- 
tes, which have simple transformation properties under 
Lorentz transformations. We start with usual multipar- 
ticle momentum states which are tensor products of ir- 
reducible representations of the Poincare group. We ob- 
serve that under a Lorentz-transformation (see Eq. ||S}) 
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U\\pi, (Jl,P2, <?2, ■ ■ ■ ,Pn, &n) = Y \Ap u a' 1 ,Ap 2 ,a' 2 ,...,Ap n ,a' n ) HD^Rwip^A)) , (12) 



where each of the D-functions depends on a different 
Wigner rotation R\y(Pi, A). This implies that one can- 
not couple angular momenta in the standard way. So it 
is desirable to have more uniform transformation proper- 
ties of n-particle states under a Lorentz transformation. 
This is the case for velocity states. For such states all spin 
projections and also the individual particle momenta are 
subject to the same Wigner rotation and the effect of 
the Lorentz transformation goes mainly into the overall 
velocity v. We now show the construction of velocity 
states in detail and make their Lorentz transformation 
properties evident. 

We consider an n-particle system with particle mo- 
menta pi and spin projections cr^, i = 1, . . . , n, and start 
by defining internal momenta ki via 

h = B- l {v) Pl , (13) 



This equation makes evident that a velocity state is a lin- 
ear combination of multiparticle momentum states. We 
also note that velocity states transform irreducibly un- 
der transformations of the Poincare group. Concerning 
notation we in general write Ui to denote spin projection 
variables, but for velocity states and when using internal 
variables of a system we will instead write /ii for spin 
projections to make a clear distinction between general 
and internal variables. In this sense the <7j should al- 
ways appear together with the pi , whereas the pi appear 
together with the ki. 

Next we study the Lorentz-transformation properties 
of a velocity state. We again apply the general boost 
operator U\ to the velocity state (|15ll and combine its 
action with that of U B j v y Using Eq. @ one obtains 



I 

where B c (v) is a canonical spin boost, i. e. a rotation- 
less Lorentz transformation, which transforms our sys- 
tem from its rest frame to total velocity v. The momenta 
k\ , k 2 , . . . , k n satisfy 

n 

£*i = (14) 

i=l 

and thus only n — 1 of them are linearly independent. 
The remaining independent variable is the overall four 
velocity v of the system. 

The construction of a velocity state can be viewed as 
starting from a multiparticle momentum state in its rest 
frame. This state is then boosted to overall velocity v by 
means of the canonical spin boost whose inverse is used 
in Eq. (|13fl to yield the velocity state 



(15) 
(16) 

I 



UA\v,ki,fJ,i,k2,fi 2 , ...,&„, Mn) = 

= UfJJB c (v)\kl, Pi, k 2 ,fJ-2, ■ ■ ■ , k n , p n ) 

= U ABc (v)\kl, Vl, k 2 ,fl 2 , ■ ■ ■ , k n , Pn) 

= U Bc {kv) iJivl^li Mil ^2, P2, ■ ■ ■ , fcn,Mn) 

= U Bc(Av - ) U Rw \k 1 ,p 1 , k 2 , P2 , ■ ■ ■ , k n i A*n > (17) 
with the Wigner rotation 

Rw = B~\Av)AB c {v) . (18) 

Evaluating the action of Ur w on the velocity state, we 
get 



\v, h, pi, k 2 , P2, ■ ■ ■ , k n ,fx n ) ■— U Bc ( v )\kx, pi, k 2 ,fi 2 , ■ ■ ■ , k n , p n ) 

n ^ 

= ^ \PU<7uP2,<r2i ■ ■ ■ )Pn,°n) JJ Dj ilH (Rw(k, B c (v))) 

(7 1 ,CT2 ,...,(T n — ±-| 4 — 1 
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Ua\v, h, hi, fc 2 ,M2, ■ ■ • , k n , fi n ) = U B ^ Av ^U Rw \ki, fii,k 2 , ^2, ■ ■ ■ ,k n , fj, n ) 

n 

= t r s c (Ai)) \Rwki,H'iiRwk2,^2, ■ ■ ■ iRwk n ,V n )X\D^,^{R w ) 

n ^ 

= ^2 \^v,Rwki,H[,Rwk2,^2,---,Rwk n ,^n)Y[D^^(Rw)- (19) 



In this derivation use has been made of the fact that 
for canonical spin boosts the Wigner rotation Ryy cor- 
responding to a rotation R is the rotation itself. It is 
helpful to notice here that a rotation is also a Lorentz 
transformation. One can now clearly see that the rota- 
tion appearing in the D-functions and in the state is the 
same for all Hi an d all ki. So one can couple spins and 
also orbital angular momenta using the standard addi- 
tion rules, which was the desired goal. 



III. MASS OPERATOR 

In this section we discuss properties of the mass oper- 
ator and its role in our approach. 

In Refs. 0, Q one can find a general procedure for 
adding interactions to a system of free relativistic parti- 
cles so that Poincare invariance is preserved. Such a pro- 
cedure, called the Bakamjian-Thomas construction, adds 
an interaction to the free mass operator. The Bakamjian- 
Thomas construction in the point form involves the (free) 
four velocity operator V which is introduced by express- 
ing the free four-momentum operator as 



p'' 
o 



M V£ 



(20) 



Interactions are added by perturbing the free mass 
operator, Mo — > M = Mq + Mj, in such a way that 
Eqs. (|4I5[) are satisfied (i. e. the components of the four 
momentum must commute with each other and they have 
to transform as the components of a four vector under 
Lorentz transformations) . We emphasize once more that 
this formulation reflects the fact that interactions do not 
enter the Lorentz generators, but solely the components 
of the four momentum. The linear constraints on the 
interacting part of the mass operator are that it should 
be a Lorentz scalar and commute with the free four ve- 
locity. The interacting four momentum operator is then 
reconstructed by 



P^ = M 



(21) 



where M contains the interactions and Vq is still kine- 
matical. 

Before starting to construct explicitly an interacting 
mass operator, we will discuss the free four-momentum 
and mass operators. Let us first examine the effects of the 



free four-momentum operator Pq, h = 0,1,2,3, on the 
velocity states defined in Sec. |H] We consider a system 
of n free particles with masses m„ internal momenta ki, 
overall four velocity v and spin projections Hi- Then we 
define 



z=l 



(22) 



which is the free relativistic mass of the system. Recalling 
the action of the free four-momentum operator P^ on 
usual n-particle momentum states, using equations 115|l . 
(J5J and l|22|) . and evaluating the boost explicitly we get 

Po\v,ki, /il.fo, M2,- ■ ■ ,k n , Hn) = 

= MnV tl \v,ki,Hiik 2 ,H2,---,k n ,Hn) ■ (23) 

Hence a velocity state \v, k±, hi, k2, H2, ■ ■ ■ , k n , fj, n ) is an 
eigenstate of Pq with the eigenvalue M n v^ ■ Thus we 
can split Pq according to Eq. ifflTjl and a velocity state 
\v, ki, Hi, fcji A*2i • • • i k n , ix n ) becomes also an eigenstate of 
Mo and Vq with the eigenvalues M. n and respec- 
tively: 

M \v,ki,Hi, - ■ ■ ,k n ,Hn) = M n \ki,Hi,---,k n ,Hn) 
V£ \v,ki,m, ...,k n ,Hn) = \ki, fj-i, . . . , k n , Hn) 

(24) 

The decomposition of Eq. (1201) is the starting point of 
the Bakamjian-Thomas construction in point form. This 
result is also consistent with the common definition of 
the mass operator 



(25) 



since the square of the four velocity is always the identity. 



IV. BAKAMJIAN-THOMAS TYPE VERTEX 
INTERACTION 

In this section we will show how an interacting mass 
operator that couples i and (i + 1) particle channels can 
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be derived from a field theoretical vertex interaction such 
that it fits into the Bakamjian-Thomas framework. Our 
presentation takes up the procedure suggested in Ref. 
to which we also refer for further details. We will set up a 
model with a finite number of effective degrees of freedom 
and consider the dynamical equation that describes this 
system of interacting particles. 

On a direct-sum Hilbert space for i and i + 1 particles 
the (full interacting) mass operator M becomes a matrix 
operator 

^'=(oi(J + (o)' ™ 

The two subspaces are coupled by the vertex operator K 
and T)\ denotes a free (indicated by the superscript °) 
i-particle operator, which corresponds to the relativistic 
i-particle mass l|22l) . In order to obtain an expression 

I 



The matrix element on the right-hand side of Eq. Ij29() has 
to be understood such that two particles on the "left" and 
one on the "right" are coupled by Hi(Q) in all possible 
ways. The remaining particles yield spectator conditions. 
This means that Mj is constructed from matrix elements 
of 7i/(0) between velocity states with the same overall 
velocity which furthermore can be taken to be zero since 
Til (0) is a Lorentz scalar. What we have neglected in this 
kind of construction as compared to the full interacting 
field theory with vertex interaction Tij are off-diagonal 
terms in the overall velocity, which should not occur in 
a Bakamjian-Thomas-type mass operator. Part of such 
terms can be simulated with an appropriate choice of the 

I 



with u>i = sj mf + kf denoting the particle energies and 
Am = uj[ + lu' 2 — uj\ — ll>2 — UJ3. TLj{Q) is the pseudoscalar 
interaction Hamiltonian density 

-i9PS ^(0)75X^(0) • £(0) , (32) 

evaluated at the space-time point x = 0. tp and 4> are the 
fermion and boson fields, the Gell-Mann flavor ma- 



for the vertex operator K we consider a field theoretical 
Hamiltonian density Tli(x), which describes a vertex in- 
teraction and is a polynomial in free fields, meaning that 

U A Hi{x)U^ =Hi(Ax) . (27) 

Then Hi(0) is a Lorentz scalar, since 

UkHMU^ 1 =Hi(0) . (28) 

These properties of Tii(0) can be used to define K. 
Taking the velocity-state representation and keeping in 
mind that the whole velocity dependence of a Bakamjian- 
Thomas type mass operator in point form is merely a fac- 
tor oc voS 3 (v' — v) (see App.^J we are led to introduce 
K via the relation 



(29) 

I 

vertex form factor /(Am) which guarantees also that the 
mass operator is a well defined operator on the Hilbert 
space. For velocity state matrix elements {v 1 = v) this 
form factor can be expressed as 

flip' - Pf] = f[(M' i+1 v' - M lV f] = /[(Am) 2 ] (30) 

since v 2 = 1. 

For our subsequent application to vector mesons the 
matrix element (|2T)|) should describe the coupling of a 
pseudoscalar meson to a quark. It has the particular 
form 



(31) 



trices, and gps is the pseudoscalar quark-meson coupling 
constant. The kincmatical factor in front of the matrix 
element (k±, fXi, fe, fJ,2, k3\H.i (0)\k'±, /4, k' 2 , fx'2) has been 
chosen such that (fci, /ii, &2, H2, fel and \k{ , fi[, k' 2 , fi' 2 ) can 
be taken as usual momentum states (with k\ +k% + k$ = 
and k[ + k' 2 = 0). 

We now turn to dynamical equations. The bound-state 



(u,fci,^i,fc 2 ,/! 2 , • • ■ , k l+1 , n l+1 \K\v' ,k[, ^^'z, n' 2 , . . . oc 

v S 3 (v- v)f(Am)(v = 0, ki,ni, fc 2 ,M2, ■ ■ • , h+i, m+i\Hi(Q)\v = 0, k[, fi[, k' 2 ,n' 2 , fc,',/4) • 



(v,ki,fii,k 2 ,fJ,2,k 3 \K\v ,k x , n x ,k2, fa) = 

= v 6 3 (v - v ) /[Am] - ^ 27r \,, 7 , N o ( fc i,Mi I fe,Ai2,fc 3 |H/(0)|fci,/ii,A:2,Ai2> 

y/{u>l +L02+ UJ3) 3 [UJi + 0J 2 Y 



r 
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problem of nonrelativistic quantum mechanics is usually 
reduced to the stationary Schrodinger equation 



(33) 



In non-relativistic quantum mechanics the Hamiltonian 
H is the only generator of the Galilei group that contains 
interactions, the remaining nine generators are kinemat- 
ical, i. e. free of interactions. In relativistic quantum me- 
chanics one deals with the Poincare group instead and 
due to the different underlying Lie algebra at least three 
generators contain interaction terms. In the point form, 
as already mentioned above, interactions are contained 
in all components of the four-momentum whereas the 
Lorentz generators are kincmatical. This means that one 
has to solve the system of dynamical equations 



(34) 



instead of Eq. (|33() to obtain the simultaneous eigen- 
states of the components of the four-momentum oper- 
ator P 1 *. Within the Bakamjian-Thomas framework the 
interaction dependence of the four-momentum operator 
becomes particularly simple (cf. Eq. iJHJ) and Eq. i|34|) 
reduces to 



M |W) = m I*) 



(35) 



We note that the v dependence of the total wave function 
of the system ^(v, ki,ni, . . .) = (v, ki,fii, . . . \ty) factors 
out after projection of Eq. (|35|l onto velocity states. 

Starting from the coupled-channel mass operator of 
Eq. I|26l) and eliminating the i + 1 particle channel one ar- 
rives at an equation for only i particles, with one-particle 
exchange between any two of them. We write down the 
eigenvalue equation for the mass operator of the two- 
channel problem: 



1*0 
*i+i! 



1*0 
*i+i) 



(36) 



where and l^i+i) are states living on the i- and 
particle subspaces, respectively. T>° and 2?° +1 denote the 
free i- and i + 1-particle masses, m is the mass eigenvalue 
of the system, and K is the vertex operator coupling the 
two channels. From these two coupled equations for 
and the latter can be eliminated to yield 



(37) 



Using velocity states, one can now turn this equation 
into an integral equation, which is a generalized eigen- 
value equation, because m also appears in the propagator 
(to - £>° +1 ) _1 - The operator ift( m - V<l +1 )- l K acts as 
an optical (one-particle exchange) potential. This optical 
potential contains, in principle, also loop contributions in 
which the exchange particle is reabsorbed by the emit- 
ting particle. But since we are interested in studying 
relativistic few-body systems which describe the dynam- 
ics of effective degrees of freedom we neglect such con- 
tributions and assume that they can be absorbed in the 
(renormalized) mass of the emitting particle. 



V. THE DYNAMICAL EQUATION 

So far we have not been very specific about the system 
we want to investigate. Equation (|36|l is a general mass 
eigenvalue equation for any coupled two-channel prob- 
lem. Up to this this point one could use this equation 
for various different systems as described at the end of 
this paper. In this section we apply it to a confined sys- 
tem of a constituent quark and a constituent antiquark 
interacting via pseudoscalar meson exchange, in order to 
give a description of vector mesons. We make use of 
the velocity-state representation, so that the Bakamjian- 
Thomas properties of the interacting mass operator can 
be fully exploited. Taking only the one-meson exchange 
dynamics into account, we could immediately start with 
Eq. with K defined according to Eqs. and 

Quarks, however, have not been observed as free par- 
ticles in nature and are therefore subject to a confining 
force. For a system of a constituent quark and a con- 
stituent antiquark interacting via pseudoscalar meson ex- 
change, i in Eq. I|37|l is equal to 2 and one has 



iVl - m)|* 2 > = K\VI - m)- x K |* 2 



(38) 



The right-hand side now corresponds to an eigenvalue- 
dependent pseudoscalar meson-exchange potential. In 
order to introduce confinement in this equation in as sim- 
ple a manner as possible we modify the relativistic kinetic 
energy terms T>\ and T>® to include harmonic oscillator 
confinement. As a next step the two-particle wave func- 
tion is expanded in terms of harmonic oscillator eigen- 
functions, which will lead to a discretization of the prob- 
lem and allow us to apply straightforward techniques for 
the numerical solution of the equation after a standard 
partial wave analysis has been carried out. 

Confinement is included in Eq. I|38|l in the diagonal 
terms in such a way that the two quarks are confined in 
the two-particle channel as well as in the three-particle 
channel, whereas the third particle, the exchange bo- 
son, is free. As we will argue in the following this still 
provides the correct Lorentz-transformation properties. 
In the two-particle channel one can introduce confine- 
ment "by hand" by substituting the free two-particle 
mass operator by a confinement one, say T>2- Such a 
confinement term a priori does not need to be of har- 
monic oscillator type. The constraint that guarantees 
Lorentz invariance is that the confinement operator must 
not depend on the overall velocity and must be a rota- 
tional scalar. In the three-particle channel only the two- 
quark subsystem is confined, which is not at rest and 
therefore has to be transformed to the correct frame. 
Using velocity states as a basis, all internal momenta 
(and also the angular momenta via the corresponding 
D-functions, see Eq. (|19l) ) are rotated by the same rota- 
tion as an effect of a Lorentz transformation. Since only 
scalar products of internal momenta appear in the diag- 
onal terms, nothing changes by an overall rotation and 
also the three-particle confinement mass operator has the 
correct Lorentz-transformation properties. 



7 



To introduce confinement in the outlined manner one 
has to make the following replacements in Eq. (|38|) 

v 2 — ► v 2 and ( 39 ) 
Vl — > VI . (40) 

The operators T>2 and T>\ are explained in more de- 
tail in App. El They are essentially square roots of 
the usual Schrodinger operator for the three-dimensional 
isotropic harmonic oscillator. Appendix 1X1 contains also 
the velocity-state representation of these operators as 
well as our actual notation for the eigenvalues. 



VI. SOLUTION METHOD AND NUMERICS 

In the preceding chapter we have encountered an equa- 
tion which has the structure of an eigenvalue equation, 
but the eigenvalue also appears in the optical poten- 
tial. Therefore one cannot directly employ standard tech- 
niques for eigenvalue problems. We use the following ap- 
proach: 

• One sets the eigenvalue in the optical potential to 
some preset value m pro and treats 

*^pre *^ Pa- 
rameter. In this way the equation becomes a linear 
eigenvalue equation and can be solved using stan- 
dard techniques to obtain eigenvalues Aj(m pre ), 
,7 = 1,2,.. ., (which, of course, depend on the value 

Of TOpro). 

• For different preset values m pre the resulting spec- 
tra will be different. 

• Interpolating the spectra for different m pre leads to 
continuous functions Xj(m pie ) of the preset eigen- 
value TO prc . The graphs of these functions exhibit 
the m pro -dependence of the spectrum caused by the 
m pre -dependence of the optical potential. 

• To find the positions of the generalized eigenvalues 
Mj, j = 1, 2, . . ., one has to solve the equations 

m = Re(Xj(m)) . (41) 

We will call Eq. (|41|) the "resonance condition" . 

• We take the real part of Xj (m) , since the eigenval- 
ues are in general (at least above the thresholds for 
the production of the various exchanged mesons) 
complex numbers. This is justified, because the 
shifts and widths resulting from the hyperfine in- 
teraction are actually perturbations to the pure os- 
cillator spectra. The imaginary part of Xj(Mj) cor- 
responds to the width of the (resonant) state with 
mass Mj, i. e. 

r(M j ) = 2|Im(A,-(M J -))|. (42) 



A procedure of this kind has, e. g., also been adopted in 
Refs. 27] and H. 

Our generalized eigenvalue equation for a two-particle 
system interacting via one-particle exchange without 
confinement is basically an integral equation. After intro- 
ducing the confinement terms in the diagonal parts of the 
coupled-channel mass operator of Eq. I|2tj|) , it is natural 
to use a harmonic-oscillator basis for the quark-antiquark 
wave function. By expressing the q-q eigenstate {^2} in 
terms of the eigenstates |w, n, I, s, j, n%j) of and tak- 
ing only a finite number of basis states the eigenvalue 
equation reduces to a system of coupled algebraic 
equations from which the mass eigenvalues m and the 
expansion coefficients A n j s ; mj (see Eq. J37j) have to 
be determined. The matrix elements of the optical poten- 
tial between the harmonic oscillator eigenstates appear in 
this set of coupled equations. These are nine-dimensional 
integrals and involve sums over various quantum num- 
bers. The integrations are done using standard Monte- 
Carlo techniques. Thereby the statistical errors have al- 
ways been kept smaller than one MeV, i. e. smaller than 
about one per mille. In the course of the numerical 
calculation the Wigner rotations for the three-particle 
intermediate-state wave functions (see Eq. (|A16|l 'l have 
been neglected. This approximation seems to be jus- 
tified by two observations. On the one hand, the pair 
of three-particle wave functions which shows up in the 
completeness relation for the three-particle intermediate 
state must have similar arguments to contribute sub- 
stantially to the nine-dimensional integral. But this also 
means that the corresponding Wigner-rotations approxi- 
mately compensate each other. On the other hand, it has 
been observed in the investi gati on of baryon form factors 
within point-form dynamics |l9j that Wigner-rotation ef- 
fects are of minor importance. The numerical effort, by 
the way, would substantially increase if Wigner rotations 
were included in our calculations. 



VII. INSTANTANEOUS APPROXIMATION 

In order to study the effects of the exchange parti- 
cle in flight as compared to the standard instantaneous 
treatment of particle exchange we perform a nonrelativis- 
tic reduction of the optical potential in the point-form 
mass operator. This is done via standard techniques and 
goes along with an "instantaneous approximation" of the 
propagator in the optical potential. "Instantaneous ap- 
proximation" means that the propagator denominator 
(2?3 — m) is replaced by the energy of the exchanged 
meson \/ q 2 + ?^m cs - ^ n ^ ne non-relativistic limit the ar- 
gument of the form factor reduces to the square of the 
three-momentum of the exchanged meson. Suppressing 
the flavor part of the hyperfine interaction, one arrives 
at the well known form for pseudoscalar meson-exchange 
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potential (see, e.g., J23) 

l/ NR (fc , fej - — — — , (43) 



47r 4toiTO2 (Z 2 + ?7i 2 



Mo 



where <f is given by 



q = k' -k 



(44) 



with fc and fc' representing cm. momenta of the incoming 
and outgoing quarks, respectively. 



VIII. MODEL PARAMETERS 



Quark masses [MeV] 




Meson masses [MeV] 


m u , m d m 3 






rax Tn n m^i 


340 500 




139 


494 547 958 


Meson-quark 


coupling 




Confinement 


gl/A-K (go/g s ) 2 


A [MeV] 


K 


a [MeV] Vo [GeV 2 ] 


0.67 1.34 


566.33 


0.81 


312 -1.04115 



TABLE I: Parameters for the point form description of vector 
mesons within the chiral constituent quark model. Apart from 
a and Vo, parameters are taken from Ref. I22T 



We adopt the parameterization of the chiral con- 
stituent quark model of Ref. for our actual calcu- 
lation of the vector-meson spectrum. In the following 
the properties and parameters of this model are briefly 
reviewed. In a constituent-quark model one deals with 
constituent quarks instead of current quarks; the con- 
stituent mass is generated dynamically and is larger than 
the corresponding current quark mass. 



m u = md = 340 MeV and 



500 MeV 



turn out to be appropriate mass values for light and 
strange (constituent) quarks. These numbers can already 
be obtained approximately from simple quark-model ar- 
guments; recent lattice calculations |3(J, |31j also hint at 
these values. The constituent quarks (and antiquarks) 
are confined and interact in addition via the exchange of 
the lightest pseudoscalar mesons which are the Goldstone 
bosons associated with chiral symmetry breaking. The 
vertex describing this interaction is constructed from the 
well-known pseudoscalar interaction Hamiltonian density 
given in Eq. (|32|l . The interaction vertex involves the 
pseudoscalar coupling constant gps between constituent 
quark and exchange meson and one cutoff parameter A^, 
i = 7r, K, rj, rj , for each meson. The parameters A^ occur 
in the meson-(anti)quark- vertex form factors fi(Am) (see 
Eqs. (|2*9"|) and (|3"Ujl ) that we are using. Following Ref. [Zj^ 
we have taken 



/.(Am) = 



A? 



A? 



Ar 



(45) 



for the functional form of these vertex form factors. The 
cutoff parameters A^ are related by 



A,; = Aq + k rrii 



(46) 



with Ao = 566.33 MeV, k — 0.81, and m.; being the mass 
of the pseudoscalar meson of type i. These vertex form 
factors go to one when Am reaches zero 1 and they go to 



1 For Am = we have four-momentum conservation at the ver- 
tex with all three particles being on-mass-shell. For t-channel 



zero like 1/Am for Am — > oo, leading to an additional 
1/ Am 2 -decay of the exchange potential. In Ref. [22] such 
a kind of form factor serves to smear out the contact 
term which occurs when (|43ll is transformed to config- 
uration space. The coupling constant g$ = gps for the 
pseudoscalar octet can be derived from the N—tt coupling 
constant via the Goldberger-Treiman relation. The value 
quoted by Glozman et al. is gpg/4-7r = 0.67. Further- 
more, two different coupling constants are used for the 
pseudoscalar meson octet and singlet, respectively. The 
ratio of the singlet to octet couplings taken in Ref. [2^ 
is (go/gs) 2 = 1-34. For our calculations the charge of 
the exchange particles is irrelevant; therefore the (small) 
mass differences between differently charged particles of 
the same sort, e.g. the ^ and the ir°, are neglected. 
The values used for the pseudoscalar meson masses are 
basically the physical masses. As in Ref. [23 we take 

m v = 140 MeV, m K = 498 MeV, m v = 547 MeV, 

and rry = 958 MeV . 



Two more parameters come from the harmonic-oscil- 
lator treatment of the quark- antiquark confinement. We 
denote the eigenvalues of T>2, i. e. the square root of the 
harmonic-oscillator eigenvalues, by 



a 2 (2n + l + -) + V Q + 4m 2 , (47) 



where a is the oscillator parameter, 4m 2 contains the rest 
masses of the quark and antiquark, and Vo leads to an 
overall shift of the spectrum (for details, see App. lA"|l . 
Since confinement is introduced in Ref. |22| in a different 
and not easily comparable way, a and Vq are free param- 
eters, a is fixed in such a way that M q and M w agree 



exchange of massive particles this can, of course, only happen 
for unphysical momenta, but it is just the kinematical situation 
(also in instant form) where the influence of the vertex form fac- 
tor is supposed to vanish and the coupling is supposed to become 
point-like. 
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with the masses of the ground state and the first excited 
state of the g spectrum. Doing this we get 

a re 312 MeV . 

This is a reasonable procedure, because the difference 
(Mio — Moo) is nearly independent of the additional hy- 
pcrfine interaction. This value for the oscillator param- 
eter a is kept fixed throughout all calculations. From 
the spectrum of the full calculations including the hyper- 
fine interaction, Vo is fixed to yield the g ground state at 
770 MeV. A suitable value for V is V = -1.04115 GeV 2 . 
All parameters of the model are summarized in Tab. [I] 

We have also done calculations without vertex form 
factors. For this purpose all parameters are kept the 
same, only Vo had to be adjusted to yield the g ground 
state at 770 MeV. One gets a slightly different value, 
namely Vo = —1.04385 GeV 2 . The calculations within 
the instantaneous approximation were performed with 
the same set of parameters as the corresponding full cal- 
culations. Finally we note that the ui and <f> flavor wave 
functions used in our calculations are the ones which cor- 
respond to ideal mixing of the singlet and octet states of 
SU{3) F . 

IX. RESULTS AND DISCUSSION 

At the beginning of this section we want to empha- 
size that our primary goal is not an optimal description 
of the meson spectrum, but rather to demonstrate with 
a simple model how the multichannel formalism devel- 
oped works and how it differs from the standard instan- 
taneous treatment of particle exchange. In our calcula- 
tions we have concentrated on the lowest-lying negative- 
parity light and strange vector mesons, i. e. mesons with 
J p = 1~ (J being the total angular momentum and P 
the parity of the system). This implies that only I = 
and I — 2 states of the harmonic-oscillator basis can con- 
tribute to the q-q wave function. Whereas I is a good 
quantum number when taking only the confining inter- 
action into account, / = and I = 2 contributions start 
to mix if the hyperfine interaction is turned on. The nu- 
merical analysis, however, reveals that the I = 2 contri- 
butions have practically no effect on the absolute masses 
(less than or at most 1 MeV, which is also the upper 
limit for our numerical accuracy). Even if compared to 
the level shift caused by the hyperfine interaction the 
I — 2 contributions are negligible with the exception of 
the two excited states of the u. For these states the I = 2 
contributions amount to 11% (first) and 18% (second ex- 
cited state) of the total level shift. In all other cases the 
contributions lack significance since they are smaller than 
the required numerical accuracy. As already explained in 
Sec. the solution of the full coupled channel problem 
involves an expansion of the vector- meson wave functions 
in terms of harmonic oscillator eigenfunctions. It turns 
out that already three basis states are enough to obtain 
convergent results on the per mille level for the ground 



and the first two excited states. For the instantaneous 
approximation of the meson exchange the convergence 
properties are worse. One needs about two times as many 
basis states as in the calculation with the full optical po- 
tential to achieve the required accuracy. It should also be 
mentioned that at those places where the harmonic oscil- 
lator eigenfunctions appear in completeness relations for 
intermediate states (cf. Eqs. <|A4(1 and l|A8[l ) the upper 
limits for the main quantum number n and the orbital 
angular momentum quantum number I have been taken 
to be the same as in the expansion of the q-q wave func- 
tion. 

The spectrum of the lowest-lying vector mesons is plot- 
ted in Fig.Q] The comparison of the full calculation and 
the pure confinement result shows that the hyperfine in- 
teraction due to (dynamical) Goldstone-boson exchange 
can be considered as a perturbation. Therefore the qual- 
itative features of the vector-meson spectrum are in our 
model essentially determined by the confinement poten- 
tial. It is thus not too surprising that only the masses of 
the ground states and the first excites states are compa- 
rable to experiment, whereas the predictions for the sec- 
ond excited states lie already much too high. To obtain 
also quantitative agreement with experiment it would 
certainly be necessary to take a confinement potential 
which is more sophisticated than our simple harmonic- 
oscillator confinement. A refined confinement potential 
which is applicable in momentum-space calculations has, 
e.g. been suggested in Ref. |33j. But as we said already 
at the beginning of this section, we rather want to study 
particle exchange within a relativistic framework and the 
conclusions about the particle exchange should not de- 
pend too much on the specific choice of the additional 
confinement potential. 

The biggest level shifts caused by the hyperfine inter- 
action are detected for the lo spectrum. This observa- 
tion can already be anticipated from the fact that the 
flavor factor at the 7r-quark vertex, the pion being also 
the lightest exchange particle, has its maximum value 
for the D meson. The u> spectrum is thus also the best 
place to study the features of our treatment of particle 
exchange. The differences between the full calculation 
and the instantaneous approximation are indeed seen to 
be most prominent in this case. Whereas the usage of a 
static meson-exchange potential for the hyperfine inter- 
action leads to an unphysically large splitting of the p and 
uj ground states, their approximate degeneracy is nearly 
preserved by our dynamical treatment of the Goldstone- 
boson exchange. The lo spectrum is obviously also most 
sensitive to the choice of the meson-quark vertex form 
factor. Comparing the results for the standard parame- 
terization of the vertex-form factors (see Sec. IVIIIjl with 
the outcome for point-like coupling, i. e. the form factors 
set to one, a striking observation can be made: whereas 
the instantaneous approximation depends very strongly 
on the form factor only a mild dependence is seen for the 
full calculation. The reason for this discrepancy is the 
difference in the propagators that make up the hyperfine 



m [MeV] 
2000- 

1800— 

1600- 



1400— 



i :•: :i 
i i 
l l 
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1000— 



800- 



600— 



EXP OSC PF IA NFF NIA 



EXP OSC PF IA NFF NIA 



EXP OSC PF IA NFF NIA 



EXP OSC PF IA NFF NIA 



FIG. 1: The spectra for the lowest-lying light and strange vector mesons. The boxes in the columns labelled by "EXP" represent the experimental values with their 
uncertainties [3^ . The other columns give our numerical results for the pure confinement interaction ("OSC"), the full calculation with dynamical mesons exchange 
("PF") and the instantaneous approximation to the meson exchange ("IA"). Corresponding results with the meson-quark vertex form factor set to one are labelled 
by "NPF" and "NIA", respectively. 
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interaction. In the instantaneous approximation it is the 
(non-relativistic) meson propagator, in the full optical 
potential it is rather the propagator of the intermedi- 
ate g-g-meson state. The q-q system in the intermediate 
state is in addition subject to confinement which acts as 
a natural cutoff and damps the dependence on the vertex 
form factors. 

Our approach does not only cover recoil effects in 
particle exchange, it provides, in principle, also non- 
perturbative predictions for vector-meson decay widths. 
As soon as the mass of a vector meson excitation be- 
comes larger than the ground state energy of the con- 
finement potential plus the mass of an exchange meson 
the corresponding channel opens and the pseudoscalar 
meson can also be emitted leaving a lower lying vector 
meson. Above such a decay threshold the optical poten- 
tial and thus the eigenvalues acquire an imaginary part 
and the width for the decay of the vector-meson reso- 
nance into the open two-particle channels can be calcu- 
lated via Eq. H42|) ■ Within our simple two-channel model 
the decay modes are restricted to p-ir, oj-tt, and p-rj. 
Among the resonances in Fig. ^ there is only one promi- 
nent, the w(1420) which decays into p-ir with a measured 
width of 174 ±40 MeV. The experimental information on 
the other resonance widths for the strong decay into one 
of the above mentioned two-particle channels is rather 
poor. Only upper bounds, which are of the order of MeV 
are given. Our theoretical results are all below 1 MeV, 
i. e. below our calculational accuracy. In the outlook we 
will discuss possible improvements of our model which 
may also lead to larger decay widths. It seems, how- 
ever, unlikely that the huge decay width of the w(1420) 
can be explained within a simple two-channel approach. 
We rather expect that other mechanisms than those in- 
cluded so far, e. g. a strong final-state interaction, have 
to be taken into account. 



X. SUMMARY AND OUTLOOK 

We have presented a Poincare invariant and Lorentz 
covariant point-form approach to the dynamical treat- 
ment of particle exchange. We have worked within the 
Bakamjian-Thomas framework, which means that the in- 
variant mass operator takes over the role of the Hamilto- 
nian in non-relativistic quantum mechanics. Operators 
and wave functions have been defined with respect to a 
velocity-state basis. Velocity states are very natural and 
advantageous for treating relativistic few-body systems 
within point-form dynamics. The starting point of our 
approach to particle exchange is a two-channel problem 
in which the i and (i + 1) particle channels are coupled via 
a vertex interaction which was derived from a field theo- 
retical Hamiltonian density such that the resulting mass 
operator is of Bakamjian-Thomas type. By reducing the 
problem to a one-channel problem for the i-particle chan- 
nel we have ended up with an optical potential which 
describes the dynamics of the particle exchange. The 



corresponding eigenvalue problem, however, is non-linear 
and has to be solved by appropriate means. Since this 
framework accounts for particle production it is able to 
provide non-perturbative predictions for (partial) decay 
widths of resonances. 

As a first application of the developed formalism 
we have investigated vector mesons within the chiral 
constituent-quark model in which the hyperfine interac- 
tion between the confined quark-antiquark pair is medi- 
ated by Goldstone-boson exchange, i. e. by the exchange 
of the lightest pseudoscalar mesons. With a simple 
harmonic-oscillator confinement and a parameterization 
of the chiral constituent-quark model that has already 
been successfully applied for the description of baryon 
spectra we have found that the hyperfine interaction 
due to Goldstone-boson exchange causes only small level 
shifts. Thus it can be considered as a perturbation of the 
confinement interaction and the confinement potential es- 
sentially determines the properties of the mass spectrum. 
The comparison of the results for the full optical poten- 
tial and the standard instantaneous meson-exchange po- 
tential revealed sizable differences, in particular for the 
uj spectrum. These differences are also reflected in the 
sensitivity to the parameterization of the meson-quark 
vertex form factors. Whereas the full calculation de- 
pends only mildly on the choice of the vertex form factors, 
the instantaneous approximation is extremely sensitive 
to changes in the form factors. Since the meson-quark 
couplings and the exchange-meson masses are subject to 
physical constraints, any reasonable parametrization of 
the Goldstone boson exchange can thus be expected to 
provide similar results in the full calculation. Our pre- 
dictions for vector-meson decay widths lie below the de- 
manded numerical accuracy and thus lack significance. 

Our conclusions from the investigation of vector me- 
sons are that a proper relativistic treatment of particle 
exchange has to go beyond the standard instantaneous 
approximation and must account for the dynamical be- 
havior of the exchange particle. The predictions for the 
vector meson spectrum could be improved with a refined 
confinement interaction. For a reasonable description of 
resonance widths it may be necessary to extend the op- 
tical potential by loop contributions, i. e. contributions 
in which the emitted meson is again absorbed by the 
emitting particle. For the present calculation we have 
assumed that such contributions go as self-energy contri- 
butions into the constituent-quark masses. But this is at 
most an approximation since the (anti) quark in a loop 
is not free, but confined. Loop contributions have, e.g., 
also been seen to be important in the semirelativistic 
treatment of the nucleon-nucleon system if one reaches 
the pion- production threshold [34( . It will be worthwhile 
and necessary to investigate their role in our coupled- 
channel formalism. This formalism should also be useful 
in treating other relativistic few-body systems which in- 
teract via particle exchange. The positronium and hydro- 
gen systems are presently under investigation. They are 
well studied within instant- and front-form dynamics and 
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APPENDIX A: NORMALIZATIONS AND MATRIX ELEMENTS 

In this appendix we collect the most important definitions and formulae used in the calculation. We start with 
some definitions concerning velocity states. Consider a system with overall four velocity v consisting of n (spin 1/2) 
fermions with masses m,. Their spins and momenta are uniquely specified by their spin projections fj,i and momenta 
ki, i = 1, . . . , n in the overall rest frame of the system. We also define uji := (mf + k 2 ) 1 / 2 . Then the completeness 
relation for the n-particle velocity states reads 

i x- f d 3 v (-\ ]3 \ (Et^) 3 



( II ^ J TT» =1 2 '- \ V >kl, Vl,k2, V2, ■ ■ ■ ,kn, Vn)(v,kl, Vl,k 2 , V2, ■ ■ ■ ,kn, Vn\ = 1 ] (Al) 

the corresponding orthogonality relation is 

(v,ki,fJ,i,k 2 ,fM2, . . . ,fc„,Mn|w',fci,Aii, ■ ■ -,k' n ,fj,' n ) = (2?r) 3 " 'r 1 — «o<5 3 (w - J| <5 3 (k; - fei)JJ W< ■ ( A2 ) 

I2^=i i=1 i=1 

The representation of the n-particle free mass operator in the basis of n-particle velocity states is 
(v,kx,(ix, . . . ,k n , VnW^v' ,k'x, y'-L, . . -,k' n ,fi' n ) = 



rin r, n— 1 n n 

(2tt) 3 " ^r 1 S v oS 3 (v - if) I] <5 3 (fc, - k>) J] ^ ^ sjm^ + fcf . . AM , 



(Ei=i w 



3=1 



A state representing the confined quark-antiquark pair is labelled by the overall velocity v of the pair and the internal 
(oscillator) quantum numbers n and I, the total spin s, as well as the total angular momentum j and its projection 
rrij . The completeness relation for such states is 



1 °° °° 1 ' +s J f d 3 v M 2 

^E E E E E / ^ ~T Ivnlsjm^ivnlsjmjl = 1 , (A4) 



(2tt) 3 ^ ^ ^ ' 

y 1 1=0 n=l s=0 j— 



\l—s\ m.j=-j 



where are just the harmonic oscillator eigenvalues (see Eq. (|A14|I 1 with main quantum number n and orbital 
angular momentum quantum number The corresponding orthogonality relation is 

2 

(v'n'l's'j'm'jlvnlsjmj) = (2tt) 3 — t v 5 (v-if) S n > n Si'iS s > s 5fjS m *. m:j . (A5) 

nl 

The mass operator T>\ for the confined pair in the above basis is 

2 

(v'n'l's'j'm'AV^vnlsjmj) = (2tt) 3 — t v S 3 (v- if) 6 n , n 8 n 8 a > s S r jS m , mj M nl . (A6) 

M nl 3 

Expansion coefficients A n i S j rrlj of the q-q wave function with respect to the harmonic-oscillator basis are defined by 
(vnlsjm^) = {vnlsjm 3 \V, ^ int ) = (2tt) = — — ^ <5 3 (^ - V) A nlsjm . (A7) 

M n l 

For our problem of (negative parity) vector mesons j = 1 and rrij = —1, 0, 1 are fixed. Furthermore, parity restricts 
spin and orbital angular momentum to s — 1 and I = 0, 2, so that the coefficients A n i S j mj depend de facto only on 
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n and I. A state describing a system of a quark- antiquark cluster and a free pseudoscalar meson is labelled by the 
quantum numbers of the cluster, the overall velocity v and the relative momentum between the cluster and the third 
particle k. The completeness and orthogonality relations for such states are 



1 °° °° f d 3 v (u! 1 + uj ) 3 

>^eEEE E E / ~ rf3K IT: o^. 3 IvKnlsjm^ivKnlsjmjl = 1 (A8) 



(2tt) 6 ^^^ n ^ J w 2uj c1 2lu 3 

(v' n'n'l' s' j'm'Avunlsjmj) = (2ir) 6 2a;cl2cj3 Wq( j3^_ ^ S 3 (k- k') 8 n , n 8vi5 s , s 6j> jS m >. mj , (A9) 



respectively. The factors in the Jacobian are defined by u c \ = yjM^+K 2 and u> 3 = ^/m 2 + k 2 . The mass operator 
T>\ for the confined pair and the free third particle in this basis takes the form 

(v' n'n'l' s' j' m'^V^vKnlsjmj) = ( 2n ) 6 J~^pj3 v o8 3 (v- &) S 3 (k- k') 8n'n8vi8 s ' s 8jij8 m i. mj (ui ci + oj 3 ) . (A10) 

The states defined above can be combined to yield the wave function for the confined quark-antiquark pair. One 
has 

= (*)! *> £ J™^ ± ± CZUCT;^ Mlfal) W fa> ■ (All) 

V mt = — l m s = — s 

The functions u n i(k) are the well-known eigenfunctions of the three-dimensional isotropic harmonic oscillator. Their 
explicit form is 



1 / 2 n + l + 2 n\ ;+i fk 2 \ /fc N ' 



"-<«=> - ^1 iwrwTTw. \y) UJ e ^ ■ (A12) 

where L n 2 is a generalized Laguerre polynomial. The corresponding normalization integral is 

[u nl (k)} 2 k 2 dk = 1 . (A13) 







The Yi mi are the usual Spherical Harmonic functions. For further details see, e.g., Refs. |35j or |3jj. The eigenvalues 
M n i of D2 are the square root of harmonic-oscillator eigenvalues, i. e. 



M n i = y/8 a 2 (2n + l + 3/2) + V + 4m 2 . (A14) 

The well-known oscillator eigenvalues have been modified by adding an overall spectral shift constant Vb and an 
averaged rest-mass term in 2 to account for the different masses of the light and strange constituent quarks. The value 
of in in this expression is easily determined for g, 0, and u>, since the masses of quark and antiquark are equal in these 
cases. For the K* we adopt an averaged mass squared of the form 



(A15) 



m 2 + m\ 

m = — - ■ 

2 

For the system of the confined quark-antiquark pair and the third free particle one can write 
(v, ki, Mi, £2, [Jv,fa\vKnlsjmj) = 

= (2ir)^ v S 3 (v- v) S 3 (k- k 3 )\ 

x E E E cz; sm c^ 2 u^Y lmi (h * 

mi=-l m s =-s fj 1 ^ 2 =±I 

x D^B-Hh/mi) B c (v 12 ) Bcih/mi)} D^B- 1 {k 2 / m 2 ) B c (v 12 ) B c (k 2 /m 2 )} , (A16) 



/ 2uj c i2uj 3 


/ 2d)i2d>2 


/ 2ui 12 2ui 3 


1 (w c i + w 3 ) 3 V 


2(u)i +w 2 ) V 


(2wi2 + 2uj 3 ) 3 
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where k — B c 1 {v X2 )ki, V12 = ( » J" 12 ] , ^12 = V m i2 + ft2 > TO i2 = + ^2, w c i = \Jm 2 cl + n 2 , D.. are Wigncr D 
functions, and cDj = y^f" 



u2 + m?. 
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